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Preface 


O cean platforms have a number of advantages over ships and buoys as ocean color data collection plat- 
forms, e.g., they are common in coastal regions, stable, frequently equipped with ample power, manned 
or frequently visited, and low cost (maintained for other purposes). During the Coastal Zone Color Scanner 
(CZCS) program, optical data was collected from the Acqua Alta Oceanographic Tower (AAOT) for algorithm 
development. At the time, the influence of a large structure on radiometric measurements, particularly ocean 
reflectance observations, had not been rigorously addressed and has cast doubt on the accuracy of such mea- 
surements. This is particularly important given the 5% water-leaving radiometric accuracy goal for satellite 
ocean color measurements. 

The SeaWiFS Project has been collaborating with researchers from the Joint Research Centre (JRC) on mea- 
surement protocols and instrumentation development using the AAOT for several years, and the collaboration 
has achieved many very important results. One particularly impressive result has been the recent demonstration 
that an automated above-water reflectance system can provide very high quality data over an extended period 
of time. This accomplishment can only be exploited if an accurate accounting of platform perturbation effects is 
available. The work documented in this report summarizes the modeling and in situ data collection and analysis 
approach that was undertaken to establish the areal extent of the perturbation and the associated correction. 
As a result, an initial methodology that can be applied to other platforms is now available. Ultimately, this 
may allow for a much more extensive, and cost-effective network of coastal validation sites. 


Greenbelt, Maryland 
April 2003 


— C. R. McClain 
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Abstract 

Large offshore structures used for the deployment of optical instruments can significantly perturb the intensity 
of the light field surrounding the optical measurement point, where different portions of the visible spectrum 
are subject to different shadowing effects. These effects degrade the quality of the acquired optical data and 
can reduce the accuracy of several derived quantities, such as those obtained by applying bio-optical algorithms 
directly to the shadow-perturbed data. As a result, optical remote sensing calibration and validation studies 
can be impaired if shadowing artifacts are not fully accounted for. In this work, the general in-water shadowing 
problem is examined for a particular case study. Backward Monte Carlo (MC) radiative transfer computations — 1 
performed in a vertically stratified, horizontally inhomogeneous, and realistic ocean-atmosphere system — are 
shown to accurately simulate the shadow-induced relative percent errors affecting the radiance and irradiance 
data profiles acquired close to an oceanographic tower. Multiparameter optical data processing has provided 
adequate representation of experimental uncertainties allowing consistent comparison with simulations. The 
more detailed simulations at the subsurface depth appear to be essentially equivalent to those obtained assuming 
a simplified ocean-atmosphere system, except in highly stratified waters. MC computations performed in 
the simplified system can be assumed, therefore, to accurately simulate the optical measurements conducted 
under more complex sampling conditions (i.e., within waters presenting moderate stratification at most). A 
previously reported correction scheme, based on the simplified MC simulations, and developed for subsurface 
shadow-removal processing of in-water optical data taken close to the investigated oceanographic tower, is then 
validated adequately under most experimental conditions. It appears feasible to generalize the present tower- 
specific approach to solve other optical sensor shadowing problems pertaining to differently shaped deployment 
platforms, and also including surrounding structures and instrument casings. 


Prologue 

A primary requirement for in-water optical data mea- 
surement activities is to minimize any perturbations neg- 
atively influencing the accuracy of the observations. This 
is particularly true when the data are used as references 
or sea truth for vicarious calibration or algorithm valida- 
tion of remotely sensed ocean color data (and derived bio- 
geochemical products). The in situ optical data accuracy 
must be defined with respect to the radiometric uncer- 
tainty budget, e.g., aimed at complying with the 5% ra- 
diometric accuracy of the Sea-viewing Wide Field-of-view 
Sensor (SeaWiFS) instrument data (Hooker and McClain 
2000 ). 

Radiometric accuracy objectives can only be achieved 
if all the optically perturbing factors are minimized, or 
accounted for, so they provide an intrinsic, or corrected, 
radiometric uncertainty that cumulatively does not exceed 
the prescribed accuracy limit. Under linear perturbation 
hypotheses, each perturbation factor must be analyzed in- 
dependently to establish its magnitude and relative impor- 
tance to the total uncertainty budget. The analyses pre- 
sented here only address perturbations originating from 
optical sensor shadowing, and emphasis is placed on the 
perturbation contribution brought by the instrument de- 
ployment structure, in this case an offshore oceanographic 
tower. 

The shadowing perturbations influencing optical mea- 
surements of the in-water radiant energy field, and in- 
duced by the instrument deployment structure or by the 


instrument casing itself, can be quantified using three- 
dimensional (3-D) backward Monte Carlo (MC) simula- 
tion techniques (Gordon 1985). The measured light field 
intensity and directional properties are different from those 
existing within an unperturbed ocean-atmosphere system, 
because of the localized modification of optical propagation 
in the vicinity of the in- water instrument and the deploy- 
ment platform. The local modification is a function of the 
spectral opacity and reflectivity, plus the finite 3-D shape, 
extension, and location of the intervening shadowing struc- 
tures with respect to the radiation source (the sun) and the 
instrument detectors. 

The accuracy of results generated by bio-optical al- 
gorithms, which provide biological and physical parame- 
ters of the investigated water body as a function of the 
in situ light field, is degraded by unaccounted for shad- 
owing artifacts. These can ultimately impair the quality 
of ocean color calibration and validation activities. Such 
negative effects must be minimized wherever possible. For 
the instrument self-shading problem, optical data correc- 
tion schemes can be based on analytical methods, as for- 
mulated by Gordon and Ding (1992) and experimentally 
validated by Zibordi and Ferrari (1995). When complex 
shadowing by large and irregular deployment structures 
must be addressed, a viable approach is customized look- 
up tables or comprehensively detailed and quasi-real time 
MC simulations, as proposed by Doyle and Zibordi (2002) 
for a specific oceanographic tower case. 
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In an initial study, Zibordi et al. (1999) considered 
three aspects of the platform perturbation problem: 

• The empirical quantification of tower-shading ef- 
fects through a first set of shadowing experiments, 

• An initial evaluation of the MC simulation capabil- 
ities and accuracy, and 

• An assessment of the correction scheme guidelines 
and development effort. 

Results from the initial inquiry showed that under ordinary 
conditions, tower-shading can reduce the absolute in- water 
radiometric values by about 1-10%, thus substantially ex- 
ceeding the allowed uncertainty in radiometric data used 
for SeaWiFS calibration and validation activities. More 
recent work (Doyle and Zibordi 2002) illustrated the ap- 
plied backward MC simulation techniques, the theoretical 
3-D validation exercises, the oceanographic tower site and 
atmosphere-ocean system modeling characteristics, as well 
as the principles and theoretical results achieved with the 
look-up table shadowing-correction scheme. 

In the investigation presented here, the in situ vali- 
dation of the proposed correction scheme is attempted by 
elaborating on the findings from a second, more dedicated, 
and extensive field campaign. The shadowing experiments 
were carefully planned to conduct accurate and spatially 
controlled radiance and irradiance measurements of the in- 
water light field existing in close proximity of the Acqua 
Alta Oceanographic Tower (AAOT), located in the north- 
ern Adriatic Sea and used regularly as an ocean color cal- 
ibration and validation site as part of the Coastal Atmos- 
phere and Sea Time-Series (CoASTS) project. The sec- 
ond campaign took place from 5-17 July 1998, with the 
sequential day of the year (SDY) covering days 186-198. 
The science team members for investigating in-water tower 
perturbations are presented in Appendix A. 

1. In Situ Sampling Equipment and Methods 

The in situ sampling equipment used during the tower- 
shading campaigns was a combination of the instruments 
normally used in the ocean color activities performed at 
the tower and those needed for the specialized experiments 
associated with quantifying in-water tower perturbations. 
The former includes a large diversity of marine and at- 
mospheric measurements for the calibration and validation 
of ocean color remote sensors, while the latter includes a 
new in-water optical system with a specialized position- 
ing capability. Tower deployments have also been used 
as an opportunity to evaluate new instruments designed 
for the special circumstances associated with the coastal 
environment. Within this objective, the tower-shading ex- 
periments were used to begin a preliminary evaluation of 
a new in-water profiler. 

2. The Long Distance Deployment System (LDDS) 

The LDDS was designed to facilitate the deployment of 
free-falling optical profilers relatively far from the AAOT. 


It consists primarily of a black plastic-covered stainless 
steel cable with one end anchored to the sea floor and the 
other end secured to the top-most deck of the AAOT. The 
LDDS can be operated in a multiple-distance or single- 
distance configuration. The multiple-distance configura- 
tion ensures accurate and quick deployment of the free- 
falling profiler from 1-35 m from the AAOT in 1 m incre- 
mental steps. This configuration was used to investigate 
the tower perturbations on in-water optical data. The 
single-distance configuration permits the deployment of a 
free-falling profiler at an arbitrary distance from the tower 
(within a maximum distance of about 35 m), although, it 
is usually used at a fixed distance of 28 m (i.e., well beyond 
the influence of perturbations associated with the tower). 
This latter configuration is used during the CoASTS mea- 
surement campaigns for the operational deployment of the 
miniature NASA Environmental Sampling System (mini- 
NESS) profiler. The semirigid frame created by the taut 
steel cable, significantly improved the capability for con- 
trolling the starting location of the deployment and the 
descent of the profiler. Departures from an ideal vertical 
cast occurred only during strong currents. 

3 . Theoretical and In Situ Data 

Analysis Methods 

The comparison of theoretical and in situ AAOT per- 
turbed radiometric data, was carried out using MC sim- 
ulations and optical profiles collected at incremental dis- 
tances from the tower superstructure. In both cases, the 
AAOT shading effects on radiometric data were quanti- 
fied with respect to far-field data (assumed unperturbed). 
The backward MC simulations of the atmosphere-ocean- 
AAOT-detector system accounted for the marine and at- 
mospheric inherent optical properties (IOPs), the seawater 
vertical stratification, the AAOT geometry, and the sen- 
sor field of view. The determination of the tower pertur- 
bations from the optical profiles required the creation of 
quality assured fixed-depth optical cross sections. These 
were then fitted with a multiparameter function to mini- 
mize noise perturbations affecting single distance measure- 
ments. Confidence bands were also produced for each fitted 
function to provide error bars suitable for the comparison 
with the theoretical data. 

4. Preliminary Results 

The measurements from three field experiments, col- 
lected under very different environmental conditions, com- 
pare well with simulated model data. Under extreme con- 
ditions, however, high uncertainties in the estimated re- 
sults are produced by excessive instrument tilt and by 
wave-focusing effects. Shadowing perturbations for very 
clear waters are restricted to within 15-20 m from the tower 
legs for upwelling radiance and downward irradiance. The 
far-field unperturbed distance is reached at 20-25 m in rel- 
atively turbid waters. Under overcast conditions, the mag- 
nitude of the shadowing effect is larger than under clear 
sky conditions and the far-field distance is reached at ap- 
proximately 30 m from the tower legs. 
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Chapter 1 


In Situ Sampling Equipment 
and Methods 

Stanford B. Hooker 
NAS A /Goddard Space Flight Center 
Greenbelt , Maryland 

Giuseppe Zibordi 

JRC/ Institute for Environment and Sustainability 
Ispra, Italy 

Abstract 

The in situ sampling equipment used during the tower-shading campaigns was a combination of the instruments 
normally used in the ocean color activities performed at the tower and those needed for the specialized experi- 
ments associated with quantifying in-water tower perturbations. The former includes a large diversity of marine 
and atmospheric measurements for the calibration and validation of ocean color remote sensors, while the latter 
includes a new in-water optical system with a specialized positioning capability. Tower deployments have also 
been used as an opportunity to evaluate new instruments designed for the special circumstances associated with 
the coastal environment. Within this objective, the tower-shading experiments were used to begin a preliminary 
evaluation of a new in-water profiler. 


1.1 INTRODUCTION 

The sampling equipment used for the second tower- 
shading campaign was a refinement of the instrumentation 
deployed during the first campaign (Zibordi et al. 1999). 
The improvements to the second field campaign were as 
follows: 

1. Use of a smaller and lighter in-water optical system 
for making profiles at a variable distance from the 
tower, and 

2. Use a more rigid and accurate deployment system 
for placing the optical system at the desired dis- 
tance from the tower. 

The new deployment system was also easier to operate, 
which allowed for more measurements to be made within 
a shorter period of time. 

The profiler that was used to measure the in-water light 
field at varying distances from the AAOT during the first 
tower-shading campaign was the Low Cost NASA Envi- 
ronmental Sampling System (LoCNESS). This profiler was 
built from off-the-shelf modular components and was de- 
signed to be a low-cost alternative to the expensive, inte- 
grated free-fall instruments at the time (Aiken et al. 1998). 
Although LoCNESS was a very capable profiler in the deep 
ocean (Hooker and Maritorena 2000), it was difficult to use 


in the shallow water in the vicinity of the AAOT. Conse- 
quently, the SeaWiFS Field Team worked with the man- 
ufacturer of LoCNESS, Satlantic, Inc. (Halifax, Canada), 
to produce a more compact version of the LoCNESS in- 
strument. 

The smaller profiler differed from the previous larger 
one by positioning both of the light sensors at the end or 
tail of the profiler, but opposite to one another on the fins. 
The significant advantage of this design is it places the two 
sensors in very nearly the same horizontal plane, which is a 
more appropriate configuration for coastal waters, wherein 
vertical complexity (usually arising from optically differing 
thin layers) is a recurring feature of the water column. This 
new layout also reduced the risk of possible sensor damage 
from an impact with the sea floor. 

For the first tower-shading campaign, the experimental 
setup began with siting a black buoy approximately 30 m 
from the southern tower leg. The buoy was aligned per- 
pendicularly to the southeastern side of the tower and dis- 
placed approximately 2 m to the side of the Wire- Stabilized 
Profiling Environmental Radiometer (WiSPER) instru- 
ment. Two pulleys were then attached — one to the buoy 
and the other to tower — and a closed loop of line (60 m 
long) with marks on it every 2.5 m was run through the 
pulleys. A cable ring was linked to one of the cable marks, 
which defined the current position of the profiler (the power 
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and telemetry cable for the profiler passed through the 
ring). The ring and, thus, the profiler, was moved to 
a selected distance from the tower leg by pulling on the 
closed loop of line, until the desired number of cable marks 
between the position of the ring and the tower leg was 
achieved. 

Although this simplistic positioning system permitted 
a reliable positioning of the optical system, it was not suf- 
ficiently rigid (and, thus, accurate) for repeated use over 
extended time periods. The continued deployment of a 
small profiler at a fixed or variable far-field distance was 
a desired capability during CoASTS campaigns, so a new 
system was designed and deployed at the AAOT. This new 
system is described in Chapt. 2. 

1.2 THE AAOT 

The AAOT is located in the northern Adriatic Sea 
(12.51°E,45.31°N) approximately 15 km southeast of the 
city of Venice. The average water depth immediately be- 
low the tower is 17.5 m and the composition of the nearby 
sea floor is primarily sand and silt. The tower was built 
in 1970 and is owned and operated by the Istituto per lo 
Studio della Dinamica delle Grandi Masse (ISDGM) of 
the Italian Consiglio Nazionale delle Ricerche\ (CNR), in 
Venice. 

The tower is composed of four levels supported by four 
large pillars. Each level is approximately 7.2 x 5.2 m 2 in size 
with the exception of the lowest level, which is 5. 2x5. 2 m 2 . 
Detailed descriptions of the AAOT and its use for optical 
field measurements are presented in Hooker et al. (1999) 
and Zibordi et al. (2002). 

1.3 IN SITU INSTRUMENTS 

The principal in situ instruments that were deployed 
to directly support the data analysis for the second tower- 
shading campaign were designed to measure the IOPs and 
apparent optical properties (AOPs) of seawater. A variety 
of other instruments, which helped characterize the biogeo- 
chemical and meteorological properties of the environment 
during optical data collection, were also used where nec- 
essary. Although the majority of the data were collected 
as continuous data sequences, either vertically or tempo- 
rally, discrete water samples were collected for laboratory 
analysis. 

1.3.1 AOP Instruments 

The miniNESS and WiSPER profilers used 7-channel 
ocean color radiance series-200 (OCR-200) sensors, as well 
as 7-channel ocean color irradiance series-200 (OCI-200) 


f The Institute for the Study of Large Masses of the Italian 
National Research Council. 


sensors. Both systems used a 16-bit analog- to-digital 
(A/D) data converter unit, the so-called DATA-100, and 
acquired data at 6 Hz. All these instruments were built by 
Satlantic, Inc. (Halifax, Canada). 

Both of the in-water systems provided profiles of up- 
welling radiance L u (z , A) and downward irradiance Ed(z , A) 
in seven spectral bands. In addition, the WiSPER Ed sen- 
sor was rotated 180° to measure the upward irradiance, 
Eu(z } X). The latter data were collected in between two 
sequential Ed profiles. 

An additional OCI-200 was used to measure the above- 
water (solar) downward irradiance Ed( 0 + ,A). The solar 
reference data were acquired simultaneously with the in- 
water profiles and were required to a) correct the in-water 
measurements for changes in the illumination conditions 
during casts, and b) derive normalized quantities, like the 
remote sensing reflectance, R T S (A). 

The center wavelengths, A, of the optical radiometers 
were closely matching the nominal center wavelengths of 
most of the SeaWiFS instrument channels: 412, 443, 490, 
510, 555, 665, and 683 nm; actual center wavelengths were 
within 1 nm of these values. A comparison of the actual 
center wavelengths for miniNESS and WiSPER is pre- 
sented in Table 1. 

Table 1. Channel numbers and center wavelengths 
(in nanometers) for the radiometric sampling sys- 
tems. The sensor systems are given with their indi- 
vidual sensor codes, which are formed from a one- 
letter designator for the type of sensor, plus a three- 
digit serial number (S/N). The solar reference, 
M099, is shown only once, but it was used with 
both profilers. All of the channels have 10 nm band- 
widths. 


Chan- 

miniNESS 

WiSPER 

nel 

R035 

1040 

M099 

R046 

1071 

i 

411.1 

411.5 

411.5 

412.3 

411.3 

2 

442.9 

442.5 

442.8 

442.8 

442.9 

3 

489.9 

489.3 

489.9 

490.5 

490.2 

4 

509.7 

509.6 

510.3 

510.8 

510.1 

5 

554.8 

555.4 

554.5 

554.9 

554.8 

6 

665.0 

665.7 

664.8 

665.8 

665.6 

7 

683.1 

683.2 

683.2 

683.9 

683.6 


A solar occulter was used to periodically shade the di- 
rect sun irradiance during Ed( 0 + , A) measurements to pro- 
vide diffuse irradiance data, Ei( 0 + ,A). The latter were 
used to determine the ratio of the diffuse to direct down- 
ward irradiance, r & , which is needed in the formulations to 
correct for instrument self-shading and AAOT perturba- 
tion effects in the underwater optical measurements. 

The absolute calibration of the miniNESS and WiS- 
PER radiometers was made just before and after the field 
measurements. An analysis of the computed calibration 
coefficients showed variations to within 1.5%. 
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1.3. 1.1 miniNESS 

The miniNESS profiler is a rocket-shaped device with 
a weighted nose and buoyant (foam) fins affixed to a tail 
bracket. The free-fall aspects of the LoCNESS and mini- 
NESS designs were derived from the SeaWiFS Free-Falling 
Advanced Light Level Sensors (SeaFALLS) profiler. The 
SeaFALLS profiler was a modified version of the Satlantic 
SeaWiFS Profiling Multichannel Radiometer (SPMR). The 
primary modification was to increase the so-called righting 
moment , which is the tendency of the nose of the profiler 
to swing downwards once the rocket is released, by adding 
larger fins and increasing the amount of weight attached 
to the nose. The larger fins have an increased surface area, 
which improves the stability of the rocket as it falls through 
the water column. 

The primary difference between miniNESS and other 
free-fall profilers is the light sensors are mounted on the 
fins (Fig. 1), rather than at the ends of the instrument 
centerline. Fin mounts allow the light sensors to be very 
nearly placed in the same horizontal plane, which is an 
important capability in waters with a complicated vertical 
structure, e.g., coastal waters. 



Fig. 1 . A schematic of the miniNESS profiler. 


The first miniNESS profiler was deployed in the deep 
ocean during a cruise from South Africa to England from 
15 May to 16 June 1998, and then again during a cruise 
from England to the Falkland Islands from 11 September 
to 17 October 1998, that is, immediately before and after 
the second tower-shading campaign. The primary reason 
for deploying miniNESS on these deep ocean cruises was to 
compare its performance against established free-fall pro- 
filers under a variety of environmental conditions. 

A comparison of miniNESS with simultaneous deploy- 
ments of LoCNESS and SeaFALLS is presented in Fig. 2. 


Although the data are well distributed with respect to the 
1:1 line, there is a noticeable bias in the 510 nm data. The 
average (over the number of simultaneous casts) unbiased 
percent difference (UPD)f between the two data sets for 
each channel ranges from —0.4 to +10.8% for the 443 and 
510 nm wavelengths, respectively. If the 510 nm channel 
is ignored, the spectrally averaged UPD value is —2.4%, 
which is to within calibration uncertainties for radiance 
sensors (Hooker et a. 2002). 



Lw( A) [|iW crrv 2 nm - 1 sr 1 ] 


Fig. 2. A comparison of water-leaving radiances 
derived from miniNESS, L^( A), versus those from 
traditional free-fall profilers, L^( A), i.e., LoCNESS 
and SeaFALLS. The dashed line is a linear fit for 
the 510 nm channel. 

Difficulties with the 510 nm channels have been doc- 
umented in other intercomparisons of Satlantic profilers 
(Hooker and Maritorena 2000), but no justification for 
the differences has ever been documented. These types 
of problems are a reminder of the continuing need to in- 
tercompare instruments and calibration procedures if un- 
certainties are to be controlled. 

1.3. 1.2 WiSPER 

WiSPER is a winched system deployed through a 
custom-built profiling rig with the optical sensors mounted 
at approximately the same depth and distance (i.e., within 
a 10 cm relative depth, and 80 cm relative distance) and 
deployed at 7.5 m from the main structure of the AAOT. 
The winch speed is approximately 0.1ms -1 , which, when 
combined with the 6 Hz data acquisition rate of the op- 
tical sensors, provides 64 measurements per meter. The 
rigidity and stability of the rig is maintained through two 


f The unbiased percent difference, tp, between two quantities 
X and Y is defined as ^ = 200(X - Y)/(X + Y). 
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taut wires, anchored between the tower and a weight on 
the sea bottom, which prevents the movement of the rig 
out of the vertical plane of the wires (Fig. 3). The wire 
stabilization and the relatively low deployment speed en- 
sures a good characterization of the subsurface water layer, 
thereby minimizing the extrapolation uncertainties in de- 
termining subsurface radiance and irradiance values. 



Fig. 3. A schematic of the WiSPER profiler. 


1.3.2 IOP Instruments 

The aforementioned complexity of the coastal water 
column can add a substantial difficulty in deriving surface 
quantities from underwater AOP measurements. The pres- 
ence of one or more optically different thin layers can al- 
ter the usual log-transformed linear decay of the measured 
near-surface optical properties, and, thus, significantly de- 
grade the accuracy of the (estimated) surface values. 

Contemporaneous IOP profiles with the AOP profiles 
provide a qualitative capability for evaluating whether or 
not the near-surface layer used for AOP extrapolations 
includes layers or intrusions that would significantly de- 
grade the ability to use a log-transformed linear fit. Con- 
sequently, the WiSPER system includes the simultaneous 
collection of spectral absorption and attenuation, a(z , A) 
and c(z, A), respectively, using a Western Environmental 
Technology Laboratories (WETLabs), Inc. (Philomath, 
Oregon) 25 cm pathlength AC-9. The nine AC-9 bands 
are in the 412-715 nm spectral range and are 10 nm wide. 
A Sea-Bird Electronics (Bellevue, Washington) 5T sub- 
mersible pump was used to provide a constant water flow 
through the AC-9. 


Zibordi et al. (1999) and Hooker et al. (2003) for previous 
tower experiments. Some methodological details are in 
order, however, because they are specific to the (in-water) 
tower-shading campaigns. 

1.4.1 AOP Methods 

The objective of the field campaign was to estimate 
the shadowing effects on in-water radiometry as a func- 
tion of distance from the AAOT structure. Briefly, each 
experiment consisted of collecting, during less than 15 min 
of time, a sequence of in-water vertical profiles of down- 
ward irradiance Ed(x, z, A, t) and upwelling nadir radiance 
L u (x, z, A,£) at varying horizontal distances x from the 
tower (with x ranging from 2.5-30 m), as a function of 
depth z, wavelength A, and acquisition time t. All optical 
data were recorded making use of OCR-200 and OCI-200 
radiometers providing radiance and irradiance readings, re- 
spectively. 

Sequential optical profiles were taken with the mini- 
NESS free-falling profiler. The miniNESS instrument was 
positioned at variable distances from the tower legs using 
the LDDS (Chapt. 2). Concurrent with miniNESS data, 
incident above-water solar irradiance data, £<*(0, 0 + , A, £), 
were collected with a Multichannel Visible Detector Sys- 
tem (MVDS). The miniNESS deployments were performed 
from the southern tower leg along the x direction, run- 
ning northwest to southeast (i.e., along a vertical sampling 
plane that intersects close to the location where the rou- 
tine WiSPER in-water optical profiles are performed, at 
7.5m from the tower legs). 

1.4.2 IOP Methods 

The IOPs used to support the analysis of the opti- 
cal data collected within the AAOT shadowing experi- 
ment were the beam attenuation c(z,A) and absorption 
a(z, A) coefficients. The coefficients c(z, A) and a(z, A) were 
computed from calibrated beam attenuation and absorp- 
tion coefficients, c t _ w (z , A) and a t _ w (z, A), respectively, ob- 
tained from the AC-9 measurements for suspended and 
dissolved optical components (not including the contribu- 
tion of pure seawater). The calibrated coefficients were 
corrected for salinity and temperature differences between 
the in situ seawater and the pure water used for labora- 
tory calibration, using the conductivity, temperature, and 
depth (CTD) profile data (WETLabs 2002). 

The measured beam attenuation coefficients corrected 
for salinity and temperature effects do not require any fur- 
ther processing, that is, 

c(z,X) = A) + c„,(A), (1) 


1.4 IN SITU METHODS 

The methods used to conduct in situ measurements 
are essentially the same as those described previously in 


where the ST superscript denotes the salinity and temper- 
ature correction, and c^A) is the beam attenuation coef- 
ficient for pure water. 
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The calibrated absorption coefficients need to be fur- 
ther corrected for scattering effects, because the finite ac- 
ceptance angle of the optics and the incomplete reflectivity 
of the absorption tube surface prevents the detector from 
collecting all the scattered light, which induces an overesti- 
mate of the retrieved absorption coefficient. In the specific 
case of the CoASTS campaigns, these perturbation effects 
were removed using the method proposed by Zaneveld et 
al. (1994). 

The Zaneveld et al. (1994) method is based on the 
removal of a variable percentage of the scattering coeffi- 
cient estimated as the difference between cfj w (z, A) and 
A). The method assumes the absorption coefficient 


of particulate and dissolved material is zero at a reference 
wavelength, A 0 = 715 nm, and the shape of the volume 
scattering function is independent of wavelength, which 
means 


a(z, A) = o" (z, A) + o„(A) - Ao), (2) 

where a w (\) is the absorption of pure water taken from 
Pope and Pry (1997) and 

y ST f .. \ \ V T (.. V X )-af- T „(z,A) ,ov 
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Chapter 2 


The Long Distance Deployment System (LDDS) 

Dirk van der Linde 

JRC /Institute for Environment and Sustainability 
Ispra, Italy 

Abstract 

The LDDS was designed to facilitate the deployment of free-falling optical profilers relatively far from the AAOT. 
It consists primarily of a black plastic-covered stainless steel cable with one end anchored to the sea floor and 
the other end secured to the top-most deck of the AAOT. The LDDS can be operated in a multiple-distance 
or single-distance configuration. The multiple-distance configuration ensures accurate and quick deployment 
of the free-falling profiler from 1-35 m from the AAOT in 1 m incremental steps. This configuration was used 
to investigate the tower perturbations on in-water optical data. The single- distance configuration permits the 
deployment of a free-falling profiler at an arbitrary distance from the tower (within a maximum distance of about 
35m), although, it is usually used at a fixed distance of 28m (i.e., well beyond the influence of perturbations 
associated with the tower). This latter configuration is used during the CoASTS measurement campaigns for 
the operational deployment of the miniNESS profiler. The semirigid frame created by the taut steel cable, 
significantly improved the capability for controlling the starting location of the deployment and the descent of 
the profiler. Departures from an ideal vertical cast occurred only during strong currents. 


2.1 INTRODUCTION 

Specialized positioning equipment was designed and 
used to deploy the miniNESS profiler for the second AAOT 
shadowing campaign. The objective was to improve the 
quality of the comparison between in situ and theoretical 
optical data, with respect to the first AAOT tower-shading 
campaign. The first tower-shading campaign (Zibordi et 
al. 1999) used a capable, but somewhat simplistic and at 
times cumbersome, cabling system to move a free-falling 
profiler to variable distances from the tower in 2.5 m incre- 
ments. With the original system, it was difficult to move 
the profiler more than 20 m away from the tower (the max- 
imum distance was 22.5 m), and one of the objectives with 
the new system was to overcome this limitation, so de- 
ployments could take place as far as 30 m away from the 
tower. 

A more accurate radiometer positioning system was 
developed to conduct in- water profiling experiments in a 
spatially more controlled and reproducible fashion. The 
final layout relied strongly on the preparatory efforts un- 
dertaken to position, as rigidly as possible, the basic com- 
ponents of the system. 

2.2 THE LDDS 

A heavy weight with a floating marker buoy was sunk 
to the sea floor approximately 80 m away from the AAOT 


southern pillar in the southeast direction. A termination 
leader with flotation was attached to the anchor before 
it was deployed. A 10 mm (black) plastic-covered, stain- 
less steel cable with flotation was attached to the float- 
ing termination leader and the other end was secured to 
a manual winch on the top deck of the tower. The steel 
cable was brought under tension by winching it tightly to 
the tower (thereby causing the cable flotation on the ter- 
mination leader and the stainless steel cable to be pulled 
underwater) . 

The taut cable extending from the anchor to the top of 
the tower constitutes the primary mechanical component 
of the LDDS design. For boating safety, the steel cable 
cannot be left installed when the platform is not in use. 
So in addition to the need for rapid installation, the steel 
cable also needed to be easily dismantled. The latter was 
accomplished by attaching flotation to the end of the steel 
cable that was secured to the termination leader. When 
the tension on the steel cable is released, the end of the 
steel cable attached to the termination leader floats to the 
surface. Once on the surface, the steel cable can be de- 
tached from the termination leader and recovered on the 
tower. The flotation on the termination leader keeps the 
termination point on the surface, so the next installation 
can proceed quickly. 

The variations in how the LDDS is used on the AAOT 
are determined by the cabling systems that are attached 


8 



J.P. Doyle, S.B. Hooker, G. Zibordi, and D. van der Linde 



83.5 m 


— .13.0 m 
- 10.0 
7.0 
-4.5 
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Fig. 4. A schematic of the experimental setup, illustrating the rigid profiling structure for WiSPER, and 
the taut- wire structure for the multiple-distance configuration of the LDDS for miniNESS. The basic system 
consists of a surface buoy (point A) attached to a heavy anchor plus a steel cable with one end attached to 
a termination leader (point B) affixed to the anchor, and the other end attached to a manual winch on the 
top deck of the tower (point C). The WiSPER system is deployed from a fixed location 7.5 m from the tower, 
and the miniNESS free-falling profiler is deployed a variable distance from the tower by moving the surfacing 
ring (point G) along a taut- wire cable extending 38 m from the tower legs as part of a triangular loop of cable 
(defined by pulleys at points D, E, and F) and operated from the work deck (point H). A cable block (point I) 
prevents the profiler from descending beyond a predetermined depth (thereby minimizing the risk of a bottom 
impact). 


to the steel cable. The multiple-distance configuration was 
designed to quickly and accurately change the deployment 
point of a free falling profiler from the AAOT up to a 
distance of 35 m in 1 m increments. The single-distance 
configuration was designed to ensure an easy deployment 
of a free- falling profiler at a single distance from the AAOT 
well clear of any perturbative effects associated with the 
tower (typically about 28 m). 

2.2.1 The Multiple-Distance Configuration 

For the multiple-distance configuration of the LDDS, 
an anchorage point and pulley is affixed at the cable-water 
intersection; this point is located approximately 38 m from 
the tower base. A second pulley is attached to the AAOT 
southern pillar about 50 cm from the water surface, and a 
third pulley is mounted vertically above the second pulley 
at a height of approximately 10 m. The latter was easily 
accessible from the third-level railing. 

A 6 mm thick black nylon rope is run through the three 
small pulleys (points D, E, and F in Fig. 4) and brought 
under tension (which creates a triangular form around the 
three pulleys). A small floating ring with an inside diame- 
ter of 30 mm is attached to the nylon rope along the base of 
the triangle (point G in Fig. 4). The ring, through which 
the instrument cable passes, acts as the deployment point. 


The nylon rope is marked in 1 m increments to accurately 
identify the deployment distance from the tower. 

An operator, standing on the end of the work deck 
(point H in Fig. 4) closest to the tower, moves the nylon 
rope (between points E and F in Fig. 4) to bring the float- 
ing ring (and the profiler) to the desired distance from the 
AAOT. When the profiler is at the desired distance, the 
nylon rope is secured and the profiler is launched by quickly 
releasing the cable held at the edge of the work deck far- 
thest from the tower. A cable block (point I in Fig. 4) that 
cannot pass through the cable ring can be used to set the 
maximum amount of cable that can be released by the op- 
erator, and, thus, the maximum depth of the cast. When 
the maximum depth is reached, the profiler is pulled to the 
surface at the point from which it was launched. 

2.2.2 The Single-Distance Configuration 

The single-distance configuration of the LDDS is shown 
in Fig. 5. It is composed of a pulley block, placed at point 
E, which includes four small guide wheels rolling over the 
main cable and a 200 mm diameter pulley guide wheel for 
the cable of the profiler. There are two small (4 mm) nylon 
lines on both sides of the pulley system. One is directly 
laced at point C; the other passes through a small pulley 
located at point D (also used with the multiple-distance 
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Fig. 5. A schematic of the single-distance configuration of the LDDS for miniNESS. The system consists 
of the same basic components used in the multiple-distance configuration (Fig. 4): a surface buoy (point A) 
attached to a heavy anchor, plus a steel cable with one end attached to a termination leader (point B) affixed 
to the anchor, and the other end attached to a manual winch on the top deck of the tower (point C). The 
unique component for the single- distance configuration is the pulley block (point E) which is used in concert 
with two nylon lines run from points E-D-C and E-C to move the pulley block (and, thus, the profiler) along 
the steel cable (points D-C). 


configuration), and is then laced at point C. When the two 
small nylon lines are used together, it is possible to position 
the pulley block, and consequently the deployment point 
E, at any location along the steel cable between points C 
and D. 

The profiler is operated from the middle deck at point 
F in Fig. 5. In the stand-by position, the profiler is kept 


about 1 m under the sea surface by blocking the power- 
telemetry cable on the railing neax point F. The maximum 
deployment depth for the profiler is marked on the cable, 
and the operator uses this mark to prevent the profiler from 
sinking too deep and impacting the bottom. The usual de- 
ployment distance, indicated by point E, is approximately 
28 m from the AAOT. 
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Chapter 3 


Theoretical and In Situ Data 
Analysis Methods 

John P. Doyle! 

JRC /Space Applications Institute 
Ispra, Italy 

Abstract 

The comparison of theoretical and in situ AAOT perturbed radiometric data, was carried out using MC sim- 
ulations and optical profiles collected at incremental distances from the tower superstructure. In both cases, 
the AAOT shading effects on radiometric data were quantified with respect to far-field data (assumed unper- 
turbed). The backward MC simulations of the atmosphere-ocean-AAOT-detector system accounted for the 
marine and atmospheric IOPs, the seawater vertical stratification, the AAOT geometry, and the sensor field of 
view. The determination of the tower perturbations from the optical profiles required the creation of quality 
assured fixed-depth optical cross sections. These were then fitted with a multiparameter function to minimize 
noise perturbations affecting single distance measurements. Confidence bands were also produced for each fitted 
function to provide error bars suitable for the comparison with the theoretical data. 


3.1 INTRODUCTION 

The comparison between simulated and experimental 
data of the AAOT perturbation effects required a) 3-D 
backward MC modeling of radiance and irradiance data 
with the specific IOPs and geometry characterizing the 
AAOT field experiments and b) the implementation of 
methods for interpreting the optical measurements focus- 
ing on the creation of fixed-depth optical cross sections 
providing an experimental quantification of the AAOT per- 
turbation effects. 

3.2 THEORETICAL METHODS 

The PHO-TRAN 3-D backward MC code, its variance 
reduction techniques, and the general atmosphere-ocean- 
AAOT-detector system are detailed in Doyle and Zibordi 
(2002). Some details of the simulation system are repeated, 
because they are novel and specific to the results presented 
here. 

3.2.1 Model Principles 

To compare in situ with simulated radiometric data, 
accurate backward MC simulations of radiative transfer 
processes in the visible and near infrared within a fully 


t Currently with the Imperial College of Science, Technology 
and Medicine, Computational Physics and Geophysics, Lon- 
don, England. 


3-D and nonhomogeneous media, were performed on dif- 
ferent atmosphere-ocean-AAOT-detector systems closely 
reproducing those characterizing the tower-shading exper- 
iments. These dedicated, realistic systems include four el- 
ements: 

1. Atmospheric and oceanic IOPs, 

2. Vertical stratification of the water column, 

3. Detailed AAOT geometry, and 

4. The actual sensor field of view together with its 
light-collecting properties. 

One simplification was the spectral reflectivity of the struc- 
ture was assumed to be zero. 

The dedicated MC simulations performed in this study 
differ from those previously described in Doyle and Zibordi 
(2002) for the inclusion of the radiometer casings in the 
modeled system. In this way, the instrument self-shading 
perturbation affecting the radiometric data were intrinsi- 
cally simulated by the MC scheme, without the need of 
applying an analytic self-shading correction (Gordon and 
Ding 1992). The geometry of the miniNESS main body 
and its attached buoyancy fins were not included in the 
simulations, because they were assumed sufficiently small 
and far from the active sensor to insignificantly perturb 
the light field at the sensor location. 

The radiometer casings are modeled as completely ab- 
sorbing cylinders (11cm high and 9cm in diameter), ver- 
tically oriented, and positioned at the required simulation 
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depth. The detectors are assumed centered on the symme- 
try axis of the radiometer casing, with point-like collecting 
areas. The field of view of the radiance sensor was 20° 
(equivalent to that of the OCT200 in water). 

3,2.2 Shading Effects 

In the following analysis, the shading effects are de- 
scribed through the percent relative error (PRE) between 
a shadow-perturbed radiometric quantity, evaluated at 
an in-water point P, and the corresponding shadow-free 
quantity, evaluated exactly as for $, but in the absence 
of the AAOT structure. The PRE is indicated by 6^, and 
is defined as 

V - 10 ° ! tt. <‘> 

where, in general, $ and may show an explicit depen- 
dence on the (x,y, z) position coordinates of ty. For the 
present system, the y-dependence was dropped, because 
the profiles are all performed along constant y. The z- 
dependence was also dropped, because the PRE values are 
evaluated at a given depth. 

Shadow-free conditions can be easily modeled through 
MC simulations; in contrast , they are difficult to determine 
experimentally, because of practical limitations associated 
with deploying a free-fall system in shallow water at the 
end of a long cable. The in situ shadow-free location is 
ideally located an infinite distance from the tower. More 
realistically, it can be located at a finite far-field distance. 
The comparison between the far-field data and those per- 
turbed by the tower at distance x can be made by assuming 
the water mass is homogeneous along its horizontal planes 
(i.e., the horizontal variation in IOPs is assumed negligi- 
ble). 

The radiometric quantities ^3 and are in-water ra- 
diances and irradiances, normalized with respect to the 
incident above- water solar irradiance data, and are used 
for the comparison of in situ versus simulated PRE results 
(i.e., for the measurement-model comparison). The nor- 
malization is introduced to account for any atmospheric 
illumination change occurring during data collection. 

3.3 IN SITU METHODS 

The experimental evaluation of the tower shading per- 
turbations relied on the use of miniNESS data deployed at 
sequential distances from the tower. The quantitative use 
of these data required three major elements: 

1. Quality assurance of the data (which included pre- 
processing and subsequent quality control proce- 
dures), 

2. Depth smoothing and binning of the radiometric 
quantities, and 

3. Application of nonlinear fitting procedures for spa- 
tial interpolation. 


The detailed aspects of these elements are presented in the 
following sections. 

The optical data were collected as a series of experi- 
ments, which were composed of a sequence of casts (Ta- 
ble 2). In most cases, the distance from the tower was 
varied, but two experiments were executed at a constant 
distance from the tower, so the contribution of environmen- 
tal variability could be assessed. The latter experiments 
were conducted at x = 7.5 m, which corresponds to the 
same distance where the WiSPER data are collected, so 
these data also allow for a good comparison with WiSPER 
data. 

Table 2. A summary of the miniNESS experi- 
ments executed in the second tower-shading cam- 
paign showing the casts involved for each exper- 
iment, the time periods covered, and the range in 
distance away from the tower, x. Experiments 7 and 
13 were executed at a constant distance of 7.5 m. 


Exp. 

Casts 

SDY 

Time [GMT] 

x [m] 

1 

2- 6 

191 

0855-0902 

3 — * 9 

2 

7- 14 

191 

1128-1135 

3-15 

3 

15- 30 

191 

1151-1220 

3-29 

4 

31- 44 

191 

1222-1235 

27 —> 3 

5 

45- 56 

191 

1245-1256 

3 — 23 

6 

57- 69 

191 

1258-1313 

23—3 

7 

70- 80 

191 

1316-1329 

7.5 

8 

81- 92 

194 

1037-1050 

3 — ► 21 

9 

93-103 

194 

1104-1112 

3-21 

10 

104-114 

194 

1113-1122 

21—3 

11 

115-125 

194 

1123-1132 

3-21 

12 

126-136 

194 

1133-1142 

21-3 

13 

137-147 

194 

1144-1153 

7.5 

14 

148-158 

194 

1216-1231 

3-21 

15 

159-169 

194 

1246-1300 

3 — 21 

16 

170-180 

194 

1301-1319 

21—3 

17 

181-191 

194 

1417-1434 

3 — 21 

18 

192-206 

195 

0807-0822 

3 — 29 

19 

207-221 

195 

0824-0846 

29—3 

20 

222-232 

195 

0946-0958 

3 — 21 


Note: Cast 1 was executed as a test of the instrument 


and the deployment system. 

3.3.1 Optical Data Processing 

The optical data processing is composed of four steps: 
a) preprocessing to create the calibrated values, b) depth 
smoothing and binning to remove noise, c) construction of 
fixed-depth profiles as a function of distance, and d) qual- 
ity control procedures. Experimental error was determined 
using nonlinear curve fitting procedures and the identifi- 
cation of confidence bands. Cast sequences, representative 
of different measurement conditions, were selected for the 
shadowing analysis and are presented in Chapt. 4, along 
with summary results. 
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3. 3. 1.1 Preprocessing 


and 


Two types of optical data were collected. The so-called 
dark or offset values were recorded with caps on the ra- 
diometers (for 3 min), and raw or signal values were col- 
lected with caps off during optical profiling. The offset 
data are denoted by the O superscript and the raw data 
by the R superscript. For example, £^(0 + ,A), E d ( 0, A), 
L° (0, A), and p°(0, A) are the offset values for the above- 
water solar irradiance, in-water downward irradiance, in- 
water upwelling radiance, and pressure, respectively. The 
offset values were averages used to establish the null or 
zero level of the sensors. In the case of the pressure sensor, 
the offset is the so-called pressure tare. 

The calibration files provided by the manufacturer (Sat- 
lantic, Inc.) for the miniNESS optical sensors included the 
spectral scale factors, denoted by the superscript 5, and 
immersion factors, denoted by the superscript I. The cal- 
ibration file also included a scale factor for the pressure 
gauge ( p s ) and coefficients for the tilt sensors. 

The net digital counts were computed as the difference 
of the signal values minus the offset values. The (linear) 
calibration coefficients were used to convert the net digi- 
tal counts into absolute (calibrated) geophysical units (Ta- 
ble 3), where the latter axe identified by the C superscript. 
The depth values of the calibrated radiometric data were 
adjusted to account for the fixed distance offset between 
the miniNESS pressure gauge and the light sensors (i.e., 
Ed and L u ). The C superscript is hereafter dropped from 
symbols referring to calibrated optical data. 


L' u (x,z,X,t) 


Lujx^ A,t) 
E d (0,0+,X,t ) 


E d ( 0, 0+,\,t o ), 


( 6 ) 


where the factor £<*(0, 0 + , A, to) is applied to preserve the 
physical units of the normalized quantities, and t 0 denotes 
the starting time for the first cast in the sequence of casts 
defining a single experiment. 


3. 3. 1.2 Smoothing and Binning 

The noise in the normalized optical profile data were 
primarily induced by wave- focusing effects and excessive 
instrument tilts. Assuming the noise is equally distributed 
between high and low values, it was reduced with a seven- 
point moving-box arithmetic averaging filter. Given the 
sampling frequency (6 Hz) and the descent speed of the 
profiler (about 0.6 ms” 1 ), each box in the filter corre- 
sponds to a vertical extent of approximately 10 cm. 

The depth coregistration of the profile data was ac- 
complished through data binning. The £^(x, z,A,£) and 
Z4(x,z, A,t) data were binned at discrete depths, z*, us- 
ing depth intervals of 0.5 m (ranging from z* — 0.25 m to 
Zi 4- 0.25 m) resulting in £^(x, z, A,£*) and L' u (x, z, A, £*) 
values, where t{ is the average time for each binning in- 
terval associated with depth z*. In-water profiles of the 
hydrographic and IOP data were binned accordingly. 


3. 3. 1.3 Fixed-Depth Profiles 


Table 3. The transformations applied to convert 
the raw (digital counts) data into absolute geophys- 
ical units. The calibrated vertical tilt of the mini- 
NESS profiler, <p c , is derived from the calibrated x 
and y two-axis tilts, and <Py, respectively. 


Symbol 

Formulation and Units 

E$ t 

Ef E* d (E? - E§) [ nWcm- 2 nm" 1 ] 

L c u 

Li Ll (Li - Li) [ M.W cm -2 nm -1 sr -1 ] 

P c 

p s ( p ft _ />) [ m ] 


Cxo + Cxi (<pi) + C x2 (<pi) 2 [°] 


CyO + + Cy 2 (<pi) 2 [°] 


tan -1 ( tan 2 ip % + tan 2 ip ^) 1 [°] 


f For the above-water solar irradiance sensor, E d ( 0 + ), 
the immersion factor is, by definition, unity. 


To remove perturbations induced by changes in the 
light field during the execution of multiple miniNESS pro- 
files, the in- water optical measurements were normalized 
with respect to the solar irradiance according to: 


E d (x,z,X,t) 


Edjpo, z, A, t) 
E d (0,0+,X,t) 


E d (0,0 + ,\,t o ), 


(5) 


The optical profiles within a single experiment were 
sequentially ordered as an increasing function of x. Op- 
tical data sequences at constant depth were then gener- 
ated from each set of the in-water profiles for each exper- 
iment. A fixed- depth cross section is hereafter referred to 
as distance-profile . 

A typical AAOT shadow- perturbed distance-profile, in 
this case taken from experiment 4 (SDY 191), is shown 
in Fig. 6a. The plot displays, for the nominal channel at 
490 nm, the E' d and the V u values as a function of x and 
at a fixed depth z = 10 m. During experiment 4 the water 
column was stratified and turbid, with little wave- focusing 
effects, and the sky was clear. The cast at x — 23 m was 
affected (at 10 m depth) by a tilt greater than 7°, so it was 
flagged and not considered. 

A general feature can be seen in Fig. 6a, which is super- 
imposed upon the underlying noise: there is a falloff in the 
radiometric signal as the sensor distance with respect to 
the tower decreases. Indeed, the AAOT structure exhibits 
a remarkable shadowing influence on the in- water radiom- 
etry: at 7.5 m distance (the deployment position for WiS- 
PER) the irradiance drops to 0.060 jxW cm” 2 nm” 1 (from 
0.070 at 27 m), which represents a 14% decrease. Over the 
same interval, the radiance value drops from 0.00060 to 
0.00045 [iW cm” 2 nm” 1 sr” 1 , which is a 25% decrease. 
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Fig. 6. Typical fixed-depth AAOT shadow-perturbed distance-profile plot for A = 490 nm as a function of 
distance x from the tower, and at a fixed depth 2 = 10 m: a) the radiometric quantities E d (dotted line and 
left scale) and L' u (solid line and right scale), and b) the fitted curves and error bars added to the radiometric 
quantities. The dashed vertical line indicates the distance, 7.5 m from the AAOT structure, where WiSPER 
is routinely deployed. 


It is important to note these shadowing estimates are 
computed using only two data points at 7.5 m, and at 27 m 
(i.e., the far- Held value). This means that any noise affect- 
ing the considered data strongly influences the computed 
shadowing perturbation effects. Similar evaluations based 
on a different far-field distance (or on the same distance 
but for a repeated experiment) could, therefore, produce 
different shadowing effects at the same investigation point. 
An improved shadowing evaluation procedure is proposed 
later on using data fitting to minimize noise effects on 
(fixed-depth) distance-profile data. 

3. 3. 1.4 Data Quality Control 

Quality control of the in situ data is performed to en- 
sure their proper selection for consistent comparison with 
modeled data. The analysis of Ed and L u depth profiles 
showed that for several experiments of the first tower- 
shading campaign (Zibordi et al. 1999), the shadowing ef- 
fects were masked by the presence of mucilage which was 
observed during many of the optical casts. 

Mucilage occurs as highly scattering and inhomoge- 
neous thick layers of gelatinous aggregates (Molin et al. 
1992), which induce spatial and temporal variability within 
the in-water light field and is a recurring feature of the 
northern Adriatic Sea in the summer months (Berthon et 
al. 1999). The multiple profiles taken in the absence of mu- 
cilage, and at a variety of distances from the tower, showed 
almost regular and correlated variations, while the multiple 


profiles taken in the presence of mucilage showed irregular 
variations with respect to both distance and depth. Ex- 
periments presenting such distance and depth irregularities 
were not included in the data analysis. 

An examination of some E d and L f u profiles showed 
higher noise for the data within the topmost water layer, 
which was probably a consequence of more pronounced 
wave effects in this layer with respect to the others. Fur- 
ther analysis of optical profiles also revealed some were de- 
graded by stratification and bottom effects, which induce 
a change as a function of depth in the normal decay of 
the optical signal. Consequently, any profiles subjected to 
these three sources of degradation were rejected for further 
analysis. 

In addition, data close to the red portion of the spec- 
trum, where the in- water signal, especially at depth, suf- 
fers more from low counts (i.e., high instrument noise) and 
from Raman scattering effects not accounted for within 
the MC simulations, were not retained. Another criteria 
applied in the quality assurance of in situ data, was the 
removal of data acquired with large tilts, i.e., <p > 7°. 

3.3.2 Experimental Uncertainties 

Experimental parameters that are important in obtain- 
ing absolute radiometric comparisons between field and 
modeled data (but were not available) are the in-water 
scattering phase function, the AAOT reflectivity, the ef- 
fective skylight distribution, and the uncertainty in sensor 
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features, e.g., the non-ideal cosine response of the irradi- 
ance collectors. These factors, however, do not generally 
vary during an experiment, so they were assumed to not 
significantly influence the intercomparison between simu- 
lated and measured data through the proposed approach 
based on the far-field normalization. 

The experimental distance-profile E d and L' u data fall 
close to smooth monotonic curves. Such ideal curves would 
be obtained experimentally if all the sources of disturbance 
(e.g., wave focusing, horizontal inhomogeneities, excessive 
sensor tilts) were removed or at least minimized. 

An error minimization procedure is, therefore, needed 
to remove the source of variability within single experi- 
ments. This procedure is illustrated in the following sec- 
tions, where a multiparameter, nonlinear curve fitting al- 
gorithm is presented as applied to the distance-profile data. 
The method applied for the determination of a nonlinear 
confidence band is also described. 

3.3. 2.1 Curve Fitting Distance-Profile Data 

The quantification of the in situ PRE at the distance 
from the tower at which the radiometric quantities can be 
assumed uncontaminated by the superstructure needs to 
be defined. The selection of one fixed reference distance, 
which is assumed unperturbed, is one way to quantify the 
experimental shadowing effects. This approach may not 
always be satisfactory, because unperturbed points are not 
always located at the same position for all experiments. It 
is, therefore, preferable to adopt a procedure, by which 
an optimal fitting curve may be drawn through all the 
experimental data, and then to evaluate PRE values from 
the fitted values. In this case, all the deployment distances 
provide equal weight to the final shadowing curve, and 
the evaluated PRE is less dependent on the noise affecting 
single data points. 

To minimize surface wave and bottom perturbations, 
the analysis of E' d and L f u distance-profiles was generally 
carried out using depth bins at intermediate depths be- 
tween the sea surface and the sea floor. 

3. 3. 2. 2 The Fitting Function 

The choice of a fitting function was made based on 
its capability of reproducing the tower shading perturba- 
tions as a function of distance. A multiparameter, non- 
linear (mechanistic) growth model, which links the inde- 
pendent variable x (position) with the dependent variable 
u) (light intensity), through the so-called monomolecular 
growth function, was chosen: 

u, = a(l - (7) 

where a, /?, and 7 are positive parameters. 

The monomolecular growth function is also called the 
Box-Lucas growth function (Draper and Smith 1981) and 
arises from the integration of 

% =->(■*-"). (*) 


where, dcu/dx is the growth rate of the light intensity sig- 
nal u) (radiance or irradiance) with respect to the position 
x. This function defines the growth rate to be directly 
proportional, through 7, to the amount of growth yet to 
be achieved, a — cj. This indicates that an asymptotic 
flattening of the curve occurs at large x values. 

The parameter a represents the upper limiting value 
to the growth function, and for x = 0, the curve starts at 
the value a(l — /3). In this model, the growth rate linearly 
decreases with increasing u ;. Because of these features, the 
growth function is assumed to be adequate to describe the 
AAOT shading effects as a function of the distance from 
the superstructure. 

3. 3. 2. 3 Fitting Procedures 

A pre-existing nonlinear, least-squares fitting routine 
(non-lin-lsq.pro), available as a part of the Interactive 
Data Language (IDL) software package from Research Sys- 
tems, Inc. (Boulder, Colorado), was adopted to iteratively 
determine the fitting parameters in (7). This routine is 
based on a Marquardt-type algorithm (Marquardt 1963), 
which uses a fitting method based on the the gradient- 
expansion algorithm (which combines the best features of 
the gradient search with the method of linearizing the fit- 
ting function). 

The iterative fitting procedure stops when the conver- 
gence of the goodness of the fit (x 2 ) is to within 0.1%. The 
results of the iterative procedure are the average best-fit 
parameter values, the standard deviations of the fitting pa^ 
rameters, and the variance-covariance matrix of the fitting 
parameters. 

The best-fit parameters are used to draw the interpo- 
lating curves (Fig. 6b). The remaining information is used, 
as described in Appendix B, to determine the error bars 
associated with the curves. Figure 6 also highlights that 
the in situ noise can be minimized with averaging. Suc- 
cessful noise reduction means a lower bias is introduced 
in the evaluation of the signal intensity as a function of 
distance, and the renormalization to a value at a fixed ref- 
erence distance is not necessary. In fact, the shadow- free 
in situ distance can be assumed to occur where the steep- 
ness of the regression curve is less than a given threshold 
value. 

Values of are given in Fig. 6b for E d and L* u as “Ed 
error” and “L u error,” respectively, and evaluated at the 
fixed 7.5 m WiSPER profiling distance. These are the PRE 
values, given by 

£*(7.5) « 4(7.5) 

_ inn ffl) ~ $(7-5) W 
V(X) 

and calculated by approximating the experimental shadow- 
free quantities with the quasi-unperturbed quantities 
computed along the extrapolated fitting curve at a distance 
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x = X at which the perturbed $ value increases by less 
than 0.1% for a 1 m distance increase, i.e., where X is such 
that e<p(X + 1) - £<p(X) < 0.1%. The £<p(7.5) in situ 
values are, therefore, dependent on X , and are obtained 
separately for each shadowing experiment, making use of 
the specific set of best-fit function parameters. 

3. 3. 2.4 Fitting Difficulties 

It is possible that high noise values at a specific dis- 
tance x, may significantly influence the determination of 
the regression curve, thereby shifting the retrieved quasi- 
unperturbed point farther from, or closer to, the tower 
than the real shadow-free point would be. To minimize 
the effects of this added source of uncertainty, a solution 
is to relax the parameters of the interpolating curves while 
keeping their values consistent. 

The problem of determining the predefined variability 
of the interpolating curve between two given points is not 
different from evaluating the error bars, to a given degree 
of confidence for all points along the fitted curve. These 
error bars are referred to as “confidence bands,” because 
the curves may be drawn within them to a given degree of 
confidence, but still respecting the fit. Evaluating the con- 
fidence bands allows for a consistent comparison with MC 
simulated data, which present inherent statistical noise. A 
confidence band approach was, therefore, undertaken. 

3. 3. 2. 5 Nonlinear Fit Confidence Bands 

Once the curve fitting algorithm has provided a con- 
verging set of parameters describing the interpolating 
curve, there remains the problem of assessing to what de- 
gree of confidence such a curve represents the actual shad- 
owing effects and, thus, identifies an error bar for the in 
situ shadowing PRE. A confidence level of 95% was cho- 
sen, so that error bars were immediately comparable with 
the 2<7 error bar on the MC statistical responses, where a 
is the estimated standard deviation of the mean radiomet- 
ric quantity from MC computations (Spanier and Gelbard 
1969). 

The basic elements describing the methods applied for 
the determination of confidence bands in the multiparam- 
eter fitting theory, are provided in Appendix C. Example 
error bars on a fitted curve are shown in Fig. 6b, where 


these error bars are defined at the 95% confidence level 
(i.e., ft = 0.05). Shadowing PREs, given in the plot for 
Ed and L u at x = 7.5 m, were calculated using the fitting 
curves. 

3.4 SUMMARY 

Tower-shading effects for the AAOT are described us- 
ing the PRE between a shadow-perturbed radiometric 
quantity, evaluated at an in-water point P, and the 
corresponding shadow- free quantity, T, evaluated exactly 
as for but in the absence of the AAOT structure. The 
comparison of simulated and in situ radiometric data is 
accomplished using backward MC simulations of radia- 
tive transfer processes of the atmosphere-ocean-AAOT- 
detector system for the former and miniNESS profile data 
for the latter. The MC simulations account for: a) at- 
mospheric and oceanic IOPs, b) vertical stratification, c) 
detailed AAOT geometry (although spectral reflectivity of 
the structure was assumed zero), and d) the actual instru- 
ment geometry and sensor field of view together with its 
light-collecting properties. 

The experimental determination of tower-shading per- 
turbations through miniNESS data relies on deployments 
at sequential distances from the tower. The quantitative 
use of these data required: a) a quality assurance of the 
observations, b) depth smoothing and binning, and c) the 
application of nonlinear fitting procedures for spatial in- 
terpolation. The optical profiles belonging to the same 
experiment, defined as a sequence of profiles at increasing 
distance from the platform, were sequentially aligned as 
a function of the distance from the main tower structure. 
Optical data sequences at constant depth were generated 
from each set of in-water profiles for each experiment. 

The fixed-depth cross sections, referred to as distance- 
profiles, showed the radiance and irradiance data closely 
followed smooth monotonic curves which increased with 
increasing distance from the tower. The curves included 
random, overlapping noise from the perturbing effects of 
wave focusing, horizontal inhomogeneities, and excessive 
sensor tilts. To minimize uncertainties induced by these 
perturbations, an error minimization procedure, based on 
a multiparameter, nonlinear curve-fitting algorithm, was 
developed for the distance-profile data. 
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Abstract 

The measurements from three field experiments, collected under very different environmental conditions, com- 
pare well with simulated model data. Under extreme conditions, however, high uncertainties in the estimated 
results are produced by excessive instrument tilt and by wave-focusing effects. Shadowing perturbations for 
very clear waters are restricted to within 15-20 m from the tower legs for upwelling radiance and downward 
irradiance. The far-field unperturbed distance is reached at 20-25 m in relatively turbid waters. Under overcast 
conditions, the magnitude of the shadowing effect is larger than under clear sky conditions and the far-field 
distance is reached at approximately 30 m from the tower legs. 


4.1 INTRODUCTION 

A selection of data taken from the experiments con- 
ducted during the second AAOT tower-shading campaign 
in July 1998 (Table 2) are used for preliminary intercom- 
parisons between field measurements and model results. 
The in situ conditions were chosen to represent, as compre- 
hensively as possible, the environmental variability found 
at the AAOT site and to investigate the extent of the influ- 
ence these diverse conditions have on the AAOT shadowing 
of radiometric data. 

The data chosen for the case studies were taken from 
three distinct experiments, each one performed on a differ- 
ent day, and selected as follows: 

1. Experiment 6 (SDY 191), stratified and highly at- 
tenuating waters, clear sky, and data parsed at 9 m 
depth for all casts; 

2. Experiment 16 (SDY 194), homogeneous and very 
clear waters, clear sky (with haze), and data parsed 
at 8 m depth for all casts; and 


t Currently with Imperial College of Science, Technology and 
Medicine, Computational Physics and Geophysics, London, 
England. 


3. Experiment 18 (SDY 195), homogeneous and at- 
tenuating waters, diffuse (overcast) sky, and data 
parsed at 2 m depth for all casts. 

The data collected during these experiments and then 
parsed at the specified fixed depths are of a quality suitable 
for intercomparison with simulations, and cover a broad 
range of environmental parameters, i.e., in- water IOPs, so- 
lar zenith and azimuth angles, atmospheric optical depths, 
and illumination conditions. Quality control on selected 
distance-profiles were carried out to discard data affected 
by excessive tilts or wave-focusing effects. 

Profiles of the in situ hydrographic (temperature and 
salinity) and IOP (spectral absorption and attenuation co- 
efficients derived from AC-9 data) variables are provided in 
Figs. 7, 8, and 9 for experiments 6, 16, and 18, respectively. 
Associated with these plots, full profiles of in situ optical 
quantities perturbed by the tower shading are provided in 
Figs. 10, 11, and 12 for the same three experiments, re- 
spectively. 

The culmination of the analysis for experiments 6, 16, 
and 18 are the fixed-depth plots of distance- profile radio- 
metric data, which are provided in Figs. 13, 14, and 15, 
respectively. All of these figures show considerable shad- 
owing effects: as x increases, i.e., as the distance from 
the tower increases, the fitted radiometric values increase 
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16 21 26 30 35 40 0 0.3 0.6 0.1 1.3 2.5 

Temperature [°C] Salinity [PSU] a [nr 1 ] c [nr 1 ] 


Fig. 7. Water column hydrographic and spectral IOP data for experiment 6 on SDY 191 (stratified and highly 
attenuating waters, plus a clear sky): a) temperature, b) salinity, c) absorption, and d) attenuation. For 
the latter two parameters, the symbols correspond to the nominal wavelengths of the measurements: 412 nm 
(circles), 490 nm (squares), 555 nm (triangles), and 650 nm (diamonds). 



16 21 26 30 35 40 0 0.3 0.6 0 1.2 2.4 

Temperature [°C] Salinity [PSU] a [nr 1 ] c [nr 1 ] 


Fig. 8. Same as Fig. 7, but for experiment 16 on SDY 194 (homogeneous and very clear waters, plus a clear 
sky with haze). 
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Fig. 11. Same as Fig. 10, but for experiment 16 on SDY 194 (homogeneous and very clear waters, plus a 
clear sky with haze). 
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Fig. 12. Same as Fig. 10, but for experiment 18 on SDY 195 (very homogeneous and attenuating waters, 
plus a diffuse sky). 
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Distance [m] Distance [m] 


Fig. 13. Spectral E' d (dashed line and left scale) and L' u (solid line and right scale) distance-profiles (each 
one similar to the plot given in Fig. 6) for experiment 6 on SDY 191 (stratified and highly attenuating waters, 
plus a clear sky) at a fixed depth of 9 m and as a function of the distance from the tower: a) 412 nm, b) 
490 nm, c) 510 nm, and d) 665 nm. Fitted curves and 95% confidence bars with shadowing PRE (and ranges) 
evaluated at 7.5 m. The vertical dashed line indicates the distance the WiSPER measurements are made. 
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Fig. 16. Same as Fig. 10, but for experiment 7 on SDY 191 (characterized by stratified and highly attenuating 
waters and clear-sky conditions). This special experiment was performed by repeating 11 different casts, all 
at a 7.5 m distance from the tower. The water layer between 5-9 m depth was subject to rapidly varying 
IOPs, and the effect of this can be seen in the large variations of the radiance data from cast to cast within 
the layer. Distance-profiles falling within such unstable depth intervals should be rejected. 


monotonically until they start to level off once the far-fielcf 
distances are reached. 

An additional experiment is presented to show the im- 
portance of water column stability and its influence on the 
reproducibility of the miniNESS measurements. Figure 16 
displays the optical profiles for this additional experiment, 
during which all casts were sequentially repeated at the 
same distance (7.5 m) from the tower legs (experiment 7 
on SDY 194). The figure shows that the data within the 5- 
9 m depth interval have a bimodal expression, in the sense 
that they follow one of two paths, whereas the data at 
other depths do not. 

The distinctive variability in the 5-9 m depth interval 
is not due to varying illumination conditions. Normaliza- 
tion is applied to each cast and, moreover, such variations 
are unlikely to occur in exact coincidence with the time 
intervals involved during the 5-9 m depth sampling. Vary- 
ing shadowing effects are also not a suitable explanation, 
because the casts are all repeated at the same location and 
in a relatively short time period. 

The most likely reason for the variations observed in 
Fig. 16 is assumed to be a change in water properties 
(produced by water shears with distinctly different IOPs), 
introducing both time and space inhomogeneities. Conse- 
quently, the data in the 5-9 m depth layer were discarded in 
all experiments performed on SDY 194 and close in time to 
experiment 7. Although this is a recurring difficulty with 


trying to execute controlled experiments in the coastal en- 
vironment, the advantages of having a stable platform (the 
tower) in an oceanic environment far outweigh any nega- 
tive consequences associated with coastal variability. 

4.2 AAOT SITE CONDITIONS 

The AAOT site is located in a frontal region that can 
be characterized by either Case-1 or Case-2f water types 
(Berthon et al. 2002). The hydrological features of the area 
are primarily influenced by the northern river discharges, 
as well as by wind and rain variability. The horizontal ho- 
mogeneity in the vicinity of the tower displays extremely 
different situations. There are time periods characterized 
by high spatial inhomogeneity and other time periods char- 
acterized by high spatial uniformity. Consequently, the in 
situ daily conditions are susceptible to very short time- 
scale variations, which are superimposed on the main sea- 
sonal evolution. 

Table 4 presents a summary of some environmental 
quantities measured at the AAOT site during the second 
tower-shading campaign. Although some of the variables 


f The optical properties of Case-1 waters are solely determined 
by the phytoplankton and its derivative products (Morel and 
Prieur 1977), whereas Case-2 optical properties are also de- 
termined by other material, e.g., from terrestrial or bottom 
origin. 
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are not immediately important to the analyses presented 
here, they are presented for completeness and to show the 
extent of certain quantities. For example, the solar zenith 
angle, 6 , spanned a wide range of 11.5-47.8°, but the wind 
speed, W, was usually low (less than 5 ms” 1 ). Note the to- 
tal chlorophyll a concentration, C a , and the concentration 
of total suspended matter, Ctsm, show the water prop- 
erties sampled included clear and turbid conditions, the 
optical consequences of which are well quantified by the 
diffuse attenuation coefficient, Kd> The sums of the par- 
ticulate and yellow substance portions of the absorption 
and beam attenuation coefficients, a p + a y and c p + c y , 
respectively. 

Table 4. A summary of some environmental quan- 
tities measured at the AAOT site during the sec- 
ond tower-shading campaign. The Kd , a p -\-a yi and 
Cp+Cy values are overall averages for 490 nm and for 
the day, where the ± values indicate the standard 
deviation in percent. The sea state is coded accord- 
ing to World Meteorological Organization (1983). 


Parameter 

SDY 191 

SDY 194 

SDY 195 

Water Type 

Case-2 

Case-1 

Case-2 

C a [mgm -3 ] 

1.6 

0.3 

0.7 

Ctsm [gm~ 3 ] 

2.8 

0.8 

1.9 

K d [m- 1 ] 

0.32 

0.08 

0.23 


± 9.4% 

± 12.5% 

± 4.3% 

a p + a y [m _1 ] 

0.16 

0.03 

0.10 


± 7.7% 

± 3.8% 

±2.0% 

Cp + Cy [m" 1 ] 

2.14 

0.55 

1.26 


± 6.7% 

± 3.1% 

±0.9% 

Cloud Cover 

0/8 

0/8f 

8/8 

Illumination 

Stable 

Stable 

Stable^ 

«[°] 

11.5-34.2 

23.2-47.8 

28.5-42.2 

Stratification 

Strong 

None 

Strong 

W [ms -1 ] 

4.8 

7.3 

2.7 

Sea State 

1 

1-2 

1 


f With haze. 1 Slightly changing. 


4.2,1 Experiment 6 

Complex water stratification is observed quite clearly 
in Fig. 7 for SDY 191, and five main discrete layers can be 
roughly identified: 0-4 m, 4-6 m, 6-8 m, 8-11 m, and 11- 
15 m. The first (top) layer is probably sediment loaded, 
due to low salinity surface waters originating from terrige- 
nous river out flow. Absorption and attenuation are both 
relatively high, which means water turbidity is generally 
high. 

In this experiment, the solar zenith angle reached al- 
most 30°, the solar azimuth angle, <f > , was 235° and the sky 
conditions were clear. The Angstrom coefficient and ex- 
ponent were 0.07 and 1.5, respectively, giving an aerosol 
optical depth, t a (A), of 0.26 and 0.19, at 412 and 510 nm, 
respectively. The atmospheric pressure was 1,016.0 hPa 
giving a (Rayleigh) molecular optical depth, t r (A), of 0.32 


and 0.16, at 412 and 510 nm, respectively. The ozone load 
was 327 Dobson units giving an ozone optical depth, r 0 (A), 
of 6.6 x 10 -4 and 1.4 x 10 -4 , at 412 and 510 nm respec- 
tively. The former optical depth values were obtained by 
applying the models and relationships referenced in Doyle 
and Zibordi (2002). 

Instrument tilts for experiment 6 are the worst for the 
three case studies, because of the presence of strong and 
varying underwater currents. The currents are also the 
cause for in-water IOP time variability, which is another 
source of in situ perturbations and inhomogeneity, as inde- 
pendently documented in experiment 7 (Fig. 16). Because 
of the time difference between IOP profiles carried out with 
the AC-9 instrument at the fixed profiling location, and 
the multiple casts performed with miniNESS, a discrep- 
ancy between assumed and actual IOPs is to be expected. 
This may become the source of measurement-model dis- 
crepancies. 

In the 9 m fixed-depth profile (Fig. 13) a relatively low 
combined tilt (about 4°) occurs at a distance of 15 m, 
while a strong peak is present in the data at all wave- 
lengths. This is probably associated with wave-focusing 
effects. This degrades the quality of the fitted curve, as 
shown by the large error bars in Fig. 13 (which are more 
pronounced for Ed than for L u ). 

In summary, clear sky conditions combined with a rel- 
atively high turbidity of the water produced radiometric 
values smoothly increasing with distance from the AAOT 
structure, and thus producing pronounced shadowing ef- 
fects. 

4.2.2 Experiment 16 

A moderate, and mainly stepwise, water stratification 
is observed in Fig. 8 for SDY 194. Three layers can be 
identified at 0-10, 10-13.5, and 13.5-15 m. The relatively 
thick homogeneous surface layer was produced by a storm 
occurring in the day preceding the experiment. Absorption 
and attenuation are both quite low, when compared to the 
average AAOT site values (Berthon et al. 2002). 

In this experiment, 9 reached about 30°, and 0 « 235°. 
Sky conditions were clear and relatively transparent: the 
Angstrom coefficient and exponent values were 0.08 and 
1.8, respectively, giving t a (A) values of 0.39 and 0.27, at 
412 and 510 nm, respectively. The latter values are about 
50% higher than for SDY 191. Atmospheric pressure was 

1016.2 hPa giving a t r (A) of 0.32 and 0.16, at 412 and 
510 nm, respectively; and the ozone load was 328 Dobson 
units giving r G ( A) values of 6.6 x 10" 4 and 1.4 x 10~ 4 , at 
412 and 510 nm, respectively (the same values as for SDY 
191). 

Instrument tilts rarely exceeded 5°, because the water 
currents around the tower were moderate. At intermediate 
distances, some spikes were observed with some data, and 
these are explained by surface gravity wave perturbations 
(the sea state was higher for experiment 16). The waves 
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have an impact on the quality of the fitted curve as can 
be seen by the relatively large error bars in Fig. 14, again 
more pronounced for Ed than for L u . 

In summary, clear sky conditions combined with a high 
water transparency (which were responsible for large wave- 
focusing perturbations), produced shadowing curves that 
saturated quickly, i,e., low shadowing effects. 

4.2.3 Experiment 18 

A stepwise IOP stratification is observed in Fig. 9 for 
SDY 195. Three water layers can be identified at 0-5 m, 
5-11 m, and ll-15m, with constant values for IOPs in the 
top layer, and smoothly decreasing and then increasing 
values in the two bottom layers. The river runoff caused by 
inland heavy rain was responsible for the sediment loaded 
top layer of freshwater characterizing data in Fig. 9. The 
previously well-mixed top layer present on SDY 194 has 
disappeared and, absorption and attenuation are now both 
relatively high in the top layer. 

Sky conditions were completely overcast, presumably 
providing an almost totally uniform (i.e., isotropic) inci- 
dent radiance distribution. The atmospheric pressure was 
1004.7 hPa and the ozone load was 306 Dobson units. The 
fixed-depth profiles in Fig. 15 are at a depth of 2 m and 
suffer from significant tilts, but they did not produce any 
degrading effects because of the uniform sky light. 

The overcast skies and relatively turbid waters pro- 
duced radiometric data very smoothly increasing with dis- 
tance from the AAOT structure. The shadowing curves, 
which saturate very slowly, produce very strong shadowing 
effects on the radiometric data. 

4.3 DATA COMPARISONS 

The MC simulations were performed within specific 
sun-atmosphere-ocean-tower-radiometer systems describ- 
ing the particular environmental conditions present in the 
field during each of the tower shading experiments. In situ 
shadowing PREs are provided along with the 95% confi- 
dence intervals determined on the nonlinear function fit- 
ting the data, while MC simulated PREs are calculated 
with the ±2cr statistical error interval (which accounts for 
the 95% probability that the average response falls within 
its bounds). 

To provide preliminary measurement-model compar- 
isons, the shadowing PRE values obtained from the simu- 
lation of radiometric quantities were compared to the cor- 
responding in situ PRE values for experiments 6, 16, and 
18. Despite the effort in producing accurate in-water pro- 
files, a small number of casts appeared as outliers, that is, 
the derived optical parameters did not follow the expected 
monotonic growth of the shadowed signal as a function of 
increasing x. Outlying data from experiment 6 at a depth 
of 9 m, were easily identified and were not accounted for in 
the data analysis. 


The Ed and L u PRE values are presented spectrally 
for each radiometric quantity by giving the fitted or es- 
timated value (ave.), along with the lowest (min.) and 
highest (max.) values computed according to the experi- 
mental confidence band interval or the simulated statistical 
spread interval. The measurement-model comparisons for 
experiment 6 are presented in Table 5 and show very good 
agreement at almost all wavelengths. In those cases where 
the agreement is degraded, it is associated with the (min. 
and max.) endpoints and not the central (ave.) value. 

Table 5. Measurement and model comparisons for 
shading experiment 6 at 9 m depth. A representa- 
tive set of wavelengths are shown rather than the 
complete set of seven measured by miniNESS. 


A and 
Results 

Type 

[%] 

L 

(A) [%] 

u 

Min. 

Ave. 

Max . 

Min. 

Ave. 

Max. 

412 Data 

8.7 

10.1 

11.6 

8.4 

12.5 

14.2 

Model 

9.6 

10.0 

10.4 

12.2 

12.7 

13.1 

490 Data 

8.1 

9.7 

11.4 

12.6 

13.8 

15.3 

Model 

8.2 

9.7 

10.7 

12.1 

13.9 

15.6 

555 Data 

5.9 

7.2 

8.7 

11.8 

12.9 

15.0 

Model 

6.8 

7.4 

7.9 

12.0 

13.1 

14.2 

665 Data 

3.9 

4.5 

5.1 

3.1 

4.0 

5.4 

Model 

4.2 

4.5 

4.7 

4.0 

4.1 

4.3 


Data from experiment 16 are optimal, because they are 
from near-ideal environmental conditions (transparent wa- 
ters and clear sky conditions) and usually minimal profiler 
tilts. A relatively high tilt, however, affected the measure- 
ments in the first cast (at 3.5 m). When this bad mea- 
surement is removed from the fit, the experimental and 
simulation data are in excellent agreement, as shown in 
Table 6. 

Table 6. Measurement and model comparisons for 
shading experiment 16 at 8 m depth, after removing 
an outlying data point from the experimental data 
processing. 


A and 
Results 
Type 

[%1 

h 

(A) [%] 

u 

Min. 

Ave. 

Max. 

Min. 

Ave. 

Max. 

412 Data 

5.7 

7.2 

8.7 

8.2 

9.9 

11.6 

Model 

7.1 

7.4 

7.7 

9.2 

9.7 

10.1 

490 Data 

3.8 

5.2 

6.6 

7.0 

8.5 

10.0 

Model 

5.2 

5.6 

6.1 

7.2 

8.3 

9.4 

555 Data 

2.8 

4.0 

5.2 

5.6 

6.9 

8.2 

Model 

4.2 

4.5 

4.8 

6.0 

6.7 

7.2 

665 Data 

1.3 

2.1 

2.9 

1.9 

2.6 

3.3 

Model 

2.3 

2.4 

2.6 

2.1 

2.5 

2.9 


Simulated data from experiment 18, the overcast situ- 
ation, suffers from the assumption of an isotropic sky ra- 
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diance distribution. Departures from the assumed isotrop- 
icity are probably responsible for the simulations signifi- 
cantly overestimating the Ed shadowing PRE. 


Table 7. Measurement and model comparisons for 
shading experiment 18 at a depth of 2 m. 


A and 
Results 
Type 

[%1 

e L W [%] 

Min . Ave. Max. 

Min. Ave . Max. 

412 Data 
Model 

10.0 11.6 13.2 

18.3 

16.7 18.6 20.6 

17.9 

490 Data 
Model 

9.7 10.8 11.9 

17.7 

19.7 21.0 22.4 

16.8 

555 Data 
Model 

9.3 10.4 11.6 

17.4 

22.4 23.9 25.3 

18.0 

665 Data 
Model 

8.3 11.2 14.0 

17.2 

16.7 19.5 22.2 

16.2 


4.4 CONCLUSIONS 

The in situ data generally compared very well with the 
simulated values, although under extreme conditions al- 
lowance must be given for in situ uncertainties, mostly due 
to instrument tilts and wave-focusing effects. Attention 
also must be given to the possible presence of a significant 
nonisotropic atmospheric radiance distribution within an 
otherwise apparently ideal overcast sky. Keeping in mind 
such deviations from the modeled conditions, the major- 
ity of results from the in situ validation performed in this 
study are generally satisfactory, and support the use of the 
developed MC code to confidently assess shadowing correc- 
tions for a broad range of environmental and illumination 


conditions typically found at the AAOT site (Berthon et 
al. 2002). 

Experimental evidence, acquired at a depth within the 
water column, indicates shadowing profiles for very clear 
waters (experiment 16) seem to flatten out for both L u and 
Ed after 15-20 m from the tower legs. This far- field plateau 
is reached at 20-25 m in more turbid waters (experiment 
6). Finally, under overcast conditions (experiment 18), the 
magnitude of the shadowing effect is significantly larger 
than under other conditions, and the far-field is reached at 
a distance of about 30 m from the AAOT legs. Following 
this study, it is advisable to deploy free-falling profilers in 
the solar direction at distances in excess of 30 m from the 
AAOT legs. Simulations confirm these findings. 

At closer distances, for example at the 7.5 m deploy- 
ment distance used for WiSPER measurements, the shad- 
owing effect is notable even under clear-sky and clear- water 
conditions: both field and simulated data at a depth of 8 m 
and a wavelength of 490 nm show that with a solar zenith 
angle of 30°, the shadowing effect is within 4% for down- 
ward irradiance and within 13% for upwelling radiance. 

Real-time Monte Carlo simulations, or look-up tables 
built from precomputed simulations, can be applied to the 
atmosphere-ocean-tower system to correct for shadowing 
perturbations in optical data collected in the proximity 
of the AAOT superstructure. Based on the validations 
presented in this study, the subsurface (0'm) simplified 
system simulations, which are those carried out within the 
shadowing correction scheme proposed by Doyle and Zi- 
bordi (2002), are expected to accurately reproduce the 
measurements at the corresponding depth (i.e., those ex- 
trapolated to the 0“m level). A good accuracy is also ex- 
pected when measurements are conducted under sampling 
conditions differing from those considered in this study. 
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Appendices 


vector of m corresponding dependent variables (the shadowed 
irradiance, for instance), q be the vector of n parameters, and 
f(u; q) y a vector function of the vector variable u and parameter 
q f i.e., the nonlinear function adopted for the fitting procedure. 
For this study, the latter is the nonlinear monomolecular growth 
function lj in (4). 

The m errors rfi made in assuming the / model to represent the 
v data are 
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C. Nonlinear Confidence Bands 

Appendix A 

In- Water Tower Perturbation Science Team 


Vi = Vi - /(u*;9i-- ••,?»), (B 1) 

and define collectively the vector error 77 for /. The sum of 
squares of these errors 77 * defines the F function, which is given 
by 

F = F(q) 


The in- water tower perturbation science team members are pre- 
sented alphabetically. 



m 


John Doyle 
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London, SW7 2AZ 

UNITED KINGDOM 
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Net : j . doyleOic .ac.uk 
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Greenbelt, Maryland 20771 

Voice: 301-286-9503 

Fax: 301-286-0268 

Net: stanCfiurdbeg.gsfc.nasa.gov 
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JRC/IES/IMW T.P. 272 
1-21020 Ispra (VA) ITALY 
Voice: 39-0-332-785-362 
Fax: 39-0-332-789-034 

Net: dirk . vanderlinde© j rc . it 

Giuseppe Zibordi 
JRC/IES/IMW T.P. 272 
1-21020 Ispra (VA) 

ITALY 

Voice: +39-0-332-785-902 
Fax: +39-0-332-789-034 

Net: giuseppe.2ib0rdiCjrc.it 

Appendix B 

Multiparameter Confidence Bands 


= 'F'b’i - . . ,q n )} 2 . 

i-1 

In terms of the fitting parameters and the formulation of the 
fitting process, the F function has to be minimized with respect 
to q to obtain a measurement-model least-squares fit. The vec- 
tor of parameters q , which make the fitting function nonlinear, 
might be bounded or constrained, but for the scope of the work 
presented in this study, they are considered to be free and in- 
dependent. 

The IDL least-squares fitting routine (non-lin-lsq.pro) was 
executed by iteratively applying a gradient-expansion algo- 
rithm. The routine returned the converged best fitted param- 
eter set q minimizing F , together with the Jacobian matrix J v 
of the errors 77 with respect to the q set evaluated at q — q. The 
Hessian matrix Hf , of the sum of squares F of the error 77 , can 
be written in a approximated form as Hf — 2 J v (where 
is the transpose of J n ), while the inverse of Hf is the so-called 
variance-covariance matrix Cf of the errors 77 with respect to 
the q set of parameters. 

Let x 2 f he the sum of squares F, evaluated at q = q and divided 
by the degrees of freedom m — n. After matrix inversion on 
the JT Jrf matrix product, thus deriving the elements of the 
variance- covariance matrix, and after evaluating such elements 
at q = < 7 , thus obtaining a matrix Cf, an unbiased estimator 
of the variance of /, at any value m = ii of the independent 
variable, is given by 


var[/(u; q)] = 2 X ) £ 




dqk 


m 


A concise summary of the basic quantities and formulations, in 
multiparameter fitting theory, required to build the confidence 
bands for the fitted curve, follows the presentation of Draper 
and Smith (1981). Only essential reasoning and operative re- 
sults are given here. 

Let 771 be the number of experimental in situ observations (here 
about 15 data pairs), and n the number of parameters to be 
considered in the fit (here 3). The number of degrees of freedom 
of the fitting procedure is, thus, m — n. Let u be the vector of 
7n independent variables (the distance from the tower), v the 


where the (^) symbol indicates determination of the partial 
derivative at q = q. Taking the square root of ( B3 ) yields 
cr f (u;q), or the standard deviation on the fitted function, de- 
termined at q = q and Ui = u. 

If q is the true parameter set, then with 100(1 — i9) confidence 
the real radiometric value at point ti is contained in a confi- 
dence interval 7/(u; q), centered around the fitted function value 
f(u; q ), and given by 

= /(«;«) ± t K cr f (u; q), (B4) 
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where $ 6 [0, 1] is an arbitrary value chosen to obtain the re- 
quired percent confidence level 100(1 — #), and t K is the 100tf/2 
percentage point of the single-tailed t-distribution with m — n 
degrees of freedom, e.g., the tabulated values in Draper and 
Smith (1981). 


Appendix C 

Nonlinear Confidence Bands 

The sum of the squares function F is a function of the fit param- 
eter elements only, and in the parameter space this function can 
be represented by the contours of a surface (Draper and Smith 
1981). If the model is linear with respect to the parameters, 
the surface contours would be ellipsoids and have a single local 
minimum, at the location defined by the least-squares estima- 
tor. If the model is nonlinear, the contours are not ellipsoids, 
but tend to be irregular, and more than one local minima may 
be found. 

Nonlinear curve fitting to experimental data tends, therefore, 
to suffer from ill-conditioning (multiple local minima may be 
found in the iterative fitting procedure, depending on the cho- 
sen initial point). The accepted way to deal with this problem 
is to consider the corresponding linearized form of the nonlin- 
ear model (the one obtained by approximating the nonlinear 
model with a Taylor expansion, which is linear with respect to 
the parameters) in the vicinity of the solution minimizing the 
sum of squares. 

Prom the linearized model, the contour providing the exact el- 
lipsoidal boundary at the required confidence level is generated, 
and then labeled as being an approximate confidence boundary 
in the nonlinear case (Draper and Smith 1981). It is important 
to note the contour determined in this fashion will be a cor- 
rect confidence contour in the nonlinear case (and will not be 
elliptical in general) , and it is only the probability level that is 
approximate. 


3-D 

A/D 

AAOT 

AOPs 

CNR 

Coasts 

CTD 

GMT 

GSFC 

IDL 

IOPs 

ISDGM 


LDDS 

LoCNESS 

MC 

miniNESS 


MVDS 

NASA 


Glossary 

Three-Dimensional 

Analog- to-Digital 

Acqua Alta Oceanographic Tower 

Apparent Optical Properties 

Consiglio Nazionale delle Ricerche 
Coastal Atmosphere and Sea Time-Series 
Conductivity, Temperature, and Depth 

Greenwich Mean Time 
Goddard Space Flight Center 

Interactive Data Language 

Inherent Optical Properties 

Istituto per lo Studio della Dinamica delle 

Grandi Masse 

Long Distance Deployment System 

Low Cost NASA Environmental Sampling 

System 

Monte Carlo 

miniature NASA Environmental Sampling 
System 

Multichannel Visible Detector System 
National Aeronautics and Space Administra- 
tion 


OCI 

OCR 

PRE 

SDY 

SeaFALLS 


SeaWiFS 

S/N 

SPMR 

UPD 

WETLabs 

WiSPER 


Ocean Color Irradiance 
Ocean Color Radiance 

Percent Relative Error 

Sequential Day of the Year 

SeaWiFS Free-Falling Advanced Light Level 

Sensor 

Sea-viewing Wide Field-of-view Sensor 
Serial Number 

Satlantic SeaWiFS Profiling Multichannel Ra- 
diometer 

Unbiased Percent Difference 

Western Environmental Technology Laborer 
tories 

Wire-Stabilized Profiling Environmental Rar 
diometer 


a(A) 

a p (A) 


a w (A) 
a y (A) 


C (A) 

C 

C a 

Cf 

c p(A) 

Ctsm 

c»(A) 


Cxi 


Cy( A) 


-ST 

c t—w 


E d ( A) 
E d ( 0 + ,A) 
E' d ( A) 
Ei( A) 
Eu( A) 

/ 

F 

H f 

I 

h 


K d ( A) 


Symbols 

Spectral absorption coefficient. 

Spectral absorption coefficient of particulate mat- 
ter. 

The AC-9 absorption coefficient for particles and 
dissolved organic matter, corrected for salinity and 
temperature effects. 

The AC-9 correction term for the absorption coeffi- 
cient. 

Spectral absorption coefficient of pure seawater. 
Spectral absorption coefficient of dissolved organic 
matter. 

Spectral beam attenuation coefficient. 

Superscript for absolute geophysical units. 
Chlorophyll a concentration. 

The covariance matrix of vector errors r?. 

Spectral beam attenuation coefficient of particulate 
matter. 

Total Suspended Matter concentration. 

Spectral beam attenuation coefficient of pure sea- 
water. 

Calibration coefficients with i — 0, . . . , 2. 
Calibration coefficients with i = 0, . . . , 2. 

Spectral beam attenuation coefficient of dissolved 
organic matter (yellow substance). 

The AC-9 beam attenuation coefficient for particles 
and dissolved organic matter, corrected for salinity 
and temperature effects. 

Spectral downward irradiance. 

Spectral above-water (solar) downward irradiance. 
Normalized downward irradiance. 

Spectral indirect (diffuse) irradiance. 

Spectral upward irradiance. 

A vector function. 

A function defined by the sum of squares of rfi . 
Hessein matrix of function F. 

Superscript for the immersion coefficients. 

The confidence interval. 

The Jacobian matrix of vector error 77 . 

The transpose of J v . 

Spectral diffuse attenuation coefficient for E d ( A). 
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L u ( A) Spectral upwelling radiance. 

T ' (A) Spectral normalized upwelling radiance. 

(A) The spectral water-leaving radiance derived from 
LoCNESS and SeaFALLS measurements. 

(A) The spectral water-leaving radiance derived from 
miniNESS measurements. 

m The number of observations. 


n The number of parameters. 

O The superscript for offset data. 
p Pressure. 

p s A scale factor for the pressure gauge. 

A shadow-free quantity. 

A far-field quasi-unperturbed quantity. 

A shadow-perturbed quantity. 

q A vector of n parameters. 

q Best fitted parameter set minimizing function F. 
q True parameter set. 

<ji The q component for observation i. 

r d The ratio of the diffuse to direct downward ir radi- 
ance. 

R r8 (A) Spectral remote sensing reflectance. 

S The superscript for scale factors. 
t Time. 

t 0 The starting time. 
ti Time for depth Zi. 

u The vector of m independent variables. 

Ui The u component for observation i. 

v The vector of m dependent variables. 

Vi v component for observation i. 

W Wind speed. 

x The distance from the AAOT tower. 

X An arbitrary x quantity. 

y Horizontal displacement. 

Y An arbitrary y quantity. 

z Depth. 

Zi Binning depth. 

a Positive parameter of the oj function. 

P Positive parameter of the u function. 

7 Positive parameter of the u function. 

Percent relative error with respect to 
£y Percent relative error with respect to 

rj The vector error for the vector function /. 

7 )i The 7] component for observation i. 

0 Solar zenith angle. 

# An arbitrary value within the 0-1 interval. 

t K 100#/2 percentage points of the t-distribution with 
m — n degrees of freedom. 

A Wavelength. 

Ao A reference wavelength. 

cr Standard deviation. 

a? The standard deviation of the fitted function. 

(A) The aerosol optical depth. 

To (A) The ozone optical depth. 

Tr(A) The (Rayleigh) molecular optical depth. 

<j> Solar azimuth angle. 


<p The vertical tilt. 

<p x The i-axis component of the vertical tilt. 

<p y The y-axis component of the vertical tilt. 

X 2 The statistic quantifying the goodness of the fit. 

X 2 p The sum of squares F. 

$ Unbiased percent difference. 
uj Monomolecular function. 
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